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Abstract 


Examination of gene functions in specific tumor types improves insight in tumorigenesis and helps design better 
treatments. Due to the rarity of histiocytic/dendritic cell sarcoma in humans, it is difficult to accrue such knowledge. 
Therefore, comparative research of these cancers in predisposed dog breeds, such as the Flatcoated retriever, can be of 
value. Histiocytic sarcoma in the dog can be grouped into a soft tissue- and visceral form. The soft tissue form at first is 
localized, while the visceral form progresses more quickly to a terminal state, which might be related to variations in gene 
expression. Microarray analyses were performed on fresh-frozen tissue from Flatcoated retrievers with either soft tissue- or 
visceral histiocytic sarcoma. Expression differences of ten most significantly differentially expressed genes were validated 
with quantitative real-time PCR (q PCR) analyses. Q PCR analyses confirmed the significantly aberrant expression of three of 
the selected genes: C6 was up-regulated; CLEC12A and CCL5 were down-regulated in the visceral histiocytic sarcoma 
compared to the soft tissue form. The findings of our study indicate that these two forms of histiocytic sarcoma in the dog 
display a variation in gene expression and warrant analysis of functional changes in the expression of those genes in these 
rare sarcomas in man. 
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be studied for the benefit of affected breeds as well as for humans 
in which the rarity of diseases such as histiocytic malignancies 
hampers scientific progress. The Flatcoated retriever (PCR) has a 
strongly increased risk for HS development [16,17]. In dogs 
[3,4,16] and humans [10,11], histiocytic malignancies are grave 
conditions though the prognosis varies between subtypes [1]. HS 


Introduction 


Histiocytic malignancies include dendritic cell- and histiocytic 
sarcomas and disseminated Langerhans cell histiocytosis [1]. In 
humans, the frequency of this group of diseases is very low [2]. 
This hampers an evidence based therapeutic approach for instance 


based on selective inhibition of specific signal transduction 
pathways. Canine histiocytic sarcoma (HS) resemble the human 
histiocytic malignancies [3-6]. As for such human cancers, canine 
HS involves the proliferation of members of both histiocytic 
lineages; dendritic cells and macrophages [5,6]. The dog is 
genetically closer related to man than mice [7,8] and the study of 
genetic changes of spontaneous cancers in the dog therefore has 
high comparative value [3,9], in particular for cancers that are 
rare in the human such as histiocytic malignancies [10,11]. 
Research could lead to more insight in the pathogenesis of this 
disease and could facilitate the identification of therapeutic targets 
valuable for both species [4,7,12—15]. 

In the dog, the limited genetic flow within breeds is responsible 
for specific breed traits including disease predispositions. However, 
the disease predispositions resulting from selective inbreeding can 
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has two common subtypes. Localized, soft tissue histiocytic 
sarcoma (STHS), manifests itself as a tumor arising in the deeply 
seated soft tissues of limbs, often in association with joints. In this 
form, chemotherapy as an adjunct to tumor resection can improve 
survival in some dogs that suffer from STHS [18]. The prognosis 
for the second form; visceral histiocytic sarcoma (WHS), which is a 
multifocal and disseminated form that is manifested in internal 
organs, uniformly is very poor [3,16,17,19,20]. There is presently 
no immunohistochemical method available to definitively differ- 
entiate between STHS and VHS and some have stated that VHS 
and STHS represent two different stages along a continuum of the 
same disease [5,21]. 

Our previous research of HS in Flatcoated retrievers has 
provided evidence that the expression of nine common genes is 
altered when comparing both STHS and VHS with normal 
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spleen; indicating a common ground for the general development 
of HS [15]. As a next step we compared the two forms of HS with 
one another. In this additional study we provide evidence that, in 
addition to common changes in STHS and VHS compared to 
healthy tissues described previously, there are some marked 
genetic differences between these two common forms of histiocytic 
malignancies. 


Materials and Methods 


The experimental protocol (ID 2007.III.08.110) was peer- 
reviewed by the scientific committee of the Department of Animals 
in Science & Society, Utrecht University, The Netherlands, and 
approved by the Animal Experiments Committee of the Academic 
Biomedical Centre, Utrecht, The Netherlands. The Animal 
Experiments Committee based its decision on ‘De Wet op de 
Dierproeven’ (The Dutch ‘Experiments on Animals Act’, 1996) 
and on the ‘Dierproevenbesluit’ (the Dutch ‘animal experiments 
decree’, 1996). Both documents are available online at http:// 
wetten.overheid.nl. 


Case recruitment and histopathological evaluation 

As mentioned in the earlier study [15], all tumor samples were 
confirmed as being spontaneously occurring histiocytic malignan- 
cies obtained from (previously untreated) Dutch family-owned 
FCR. All samples were obtained with informed owner consent. 
Case selection and obtainment of tissues were similar as mentioned 
in the previous study. As a common reference pool a multitude of 
canine organs (testis, liver, spleen, prostate, duodenum, lung, 
kidney and brain) from healthy crossbreeds (n= 8) was used [15]. 

Patient details are listed in Table 1. 


RNA isolation 

RNA isolation was performed as described previously [15]. In 
short, approximately 30 mg of frozen tumor was disrupted/lysed 
and homogenized and total RNA was isolated and treated with 


Table 1. Patient details. 
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DNase using the RNeasy mini kit (Qiagen, The Netherlands) 
according to the manufacturer’s protocol. Quantity and integrity 
were assessed with the Bioanalyzer Agilent BioAnalyzer-2100 
(Bioanalyzer, Agilent Technologies, Santa Clara, CA) in combi- 
nation with an RNA 6000 Pico-LabChip. The average RNA 
integrity number 8.5 (range: 7.2—9.8) was found to be appropriate 
[22]. RNA concentration was quantified using a NanoDrop ND- 
1000 (Isogen Life Science) spectrophotometer. 


Expression profiling 

Expression profiling was performed as described previously 
[15]. In short, RNA was labeled twice and hybridized against the 
common reference RNA on dual channel arrays, with RNA 
amplifications and labeling being performed on an automated 
system (Caliper Life Sciences NV/SA, Belgium) as described 
previously [23]. Dye swap of Cy3 and Cy5 was performed to 
reduce dye bias. Hybridizations were done on a HS4800PRO 
system supplemented with QuadChambers (Tecan Benelux 
B.V.B.A.) using 500-1000 ng labeled cRNA per channel as 
described [24]. After automated data extraction using Imagene 8.0 
(BioDiscovery), printtip Loess normalization was performed [25] 
on mean spot-intensities. Dye-bias was corrected based on a 
within-set estimate as described [26]. Data were analyzed using 
ANOVA [27] (R version 2.2.1/MAANOVA version 0.98-7) 
(http://www.r-project.org/). Briefly, both tumourgroups were 
compared through the common reference channel. P-values were 
determined by a permutation F2-test, in which residuals were 
shuffled 5000 times globally. Thus analyzed 191 gene probes with 
P<0.05 after family wise error correction (FWER) were consid- 
ered significantly changed. 


Functional annotation 

To examine whether certain pathways are over- or under- 
represented in the gene list, all genes differentially expressed 
between STHS and VHS, were introduced into DAVID (http:// 
david.abcc.nciferf.gov/). 


Name Sex Pathology AO (yrs) 
D1, TJ MN STHS 77 
D3, UH FN STHS 8 
D4, BS MN STHS 6.5 
D5, BaS M STHS 76 
D6, TV MN STHS 8.1 
D7, DV MN STHS 11 
DX, YM MN STHS 9.7 
D2, BE M VHS 9.4 
Ds, DW M VHS 9.5 
D9, JV FN VHS 8.9 
D11, BT F VHS 8.5 
D12, AG M VHS 73 
D13, TR FN VHS 79 
D14, SG F VHS 41 
DX, SC MG VHS 10 


Site(s) Microarray/PCR 
shoulder Y/Y 
shoulder Y/Y 
elbow Y/Y 
knee Y/Y 
elbow Y/Y 
shoulder Y/Y 
shoulder N/Y 
liver/spleen/Inn abd Y/Y 
spleen/Inn abd Y/N 
lung/Inn mediast Y/Y 
lung Y/Y 
lung/spleen/kidney Y/Y 
lung Y/Y 
lung/Inn mediast Y/Y 
liver/spleen N/Y 


AO: Age of onset, M: male, MN: male neutered, F: female, FN: female neutered, STHS: soft tissue (localized) histiocytic sarcoma, VHS: visceral (disseminated) histiocytic 
sarcoma, Inn abd: abdominal lymphnodes, Inn mediast: mediastinal lymph nodes. Note: For cases with VHS the site sampled for gene expression is indicated in bold 
letters N=no Y=yes 
doi:10.1371/journal.pone.0098258.t001 
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Quantitative real time PCR; gene selection 

Following the outcome of the microarray expression profiling, 
11 genes were selected; namely C-type lectin domain family 12, 
member A (CLEC/2A); C-C motif chemokine 5 Precursor (Small- 
inducible cytokine A5) (CCL5_CANFA); Asporin (4SPN); CD9 
molecule (CD9); Transketolase-like 1, transcript variant 1 (TATLJ); 
Complement component 6, transcript variant 1(C’6); $100 calcium 
binding protem A1l2 (S/00A12); Immunoglobulin J polypeptide 
(IGF); S100 calcium binding protein A8 ($/00A8); Phytanoyl-CoA 
dioxygenase, peroxisomal like (PHYH). Selection of these genes 
was based on the ones most significantly differently expressed on 
the micro-array chip and their fold changes. 

Details of the qPCR reactions and primer sequences are 
depicted in Table 2. Delta Ct method, using efficiencies between 
95.0 and 104.8%, was used for both the reference- as well as the 
target genes. 


RNA isolation and cDNA synthesis 

Besides tissues from all but one patient (of which the insufficient 
tissue remained for the validating qPCR experiment) used in the 
microarray experiment (six soft tissue HS and seven visceral HS), 
two additional samples (one soft tissue HS and one visceral HS 
that met the inclusion criteria) were added. Total RNA from these 
samples was isolated. After isolation, total RNA was treated with 
DNase using the RNeasy mini kit (Qiagen, The Netherlands) 
according to the manufacturer’s protocol. 

Reverse transcription (RT) was performed of all 20 samples in a 
80 ul reaction using 2000 ng total RNA, 16 ul iScript Reaction 
mix and 4 ul iScript Reverse Transcriptase (iScript cDNA 
Synthesis kit, Bio Rad, Veenendaal, The Netherlands). The 
mixture, contain both random hexamers and oligo-dT primers,- 
was incubated 5 min. at 25°C, 30 min. at 42°C and followed by 
5 min. at 85°C. Minus RT controls were prepared from 500 ng of 
the same RNA under the same conditions, but without addition of 
reverse transcriptase. 


Reference genes and primer development 

Nine reference genes were used as the non-regulated reference 
genes for normalization, based on their stable expression in canine 
tissue [28,29] namely ribosomal protein S19 (RPS19) hypoxanthine 
phosphoribosyltransferase (HPRT), ribosomal protein L8 (RPL8), signal 
recognition particle receptor (SRPR), and ribosomal proten L1 (RPL15), 
glucuronidase, beta (GUSB), Glyceraldehyde-3-Phosphate Dehydrogenase 
(GAPDH), Beta-2-Microglobulin (B2MG), 40S ribosomal protein S5 
(RPS5). Primers for reference genes, including their optimum 
annealing temperature are listed in Table 2. 

Using Ensembl (Ensembl 70; CanFam3.1), through annotated 
transcripts, PCR primers were designed using the Perl Primer 
software (version 2.0.0.7) according to the parameters outlined in 
the Bio-Rad i-cycler manual. The specificity of each primer pair 
was confirmed by sequencing its product and also in qPCR by 
checking the meltcurve and reaction efficiency. GeNorm [30] was 
used to establish expression stability. Amplicon sequence-reactions 
were performed using BigDye v3.1 according to the manufactur- 
er’s (Life Technologies, Bleiswijk, The Netherlands) instructions on 
an ABI3130XL and analyzed in Lasergene (version 9.1 DNAS- 
TAR) and confirmed the specificity of each amplicon. 


Quantitative PCR 

Published guidelines for the qPCR experiment were followed 
according to the MIOE guidelines [31-33]. For qPCR, the CFX 
detection system (Bio-Rad.) with SYBR green fluorophore was 
used. Reactions were performed in a total volume of 10 ul 
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containing 5 ul of 2x SYBR green super mixes (Bio-Rad 
Laboratories Ltd.), 0.5 ul of each primer at 400 nM concentra- 
tion, 0.8 ul of cDNA and 3.2 ul RNase and DNase free water as 
previously described [15,34,35]. Expression analysis was per- 
formed on sample duplicates in duplicate. A minus RT sample and 
a no template control were performed as control. Expression levels 
were based on Ct values normalized using the mean of seven of the 
nine reference genes. GAPDH and B2MG did not behave fully 
stable according to GeNorm [28]. 

A Wilcoxon rank sum test was performed to determine the 
significance of differential gene expression. All results were 
Bonferroni corrected. 


Results and Discussion 


The Microarray enabled analysis of the expression of 42,034 
features; however, since only 21,682 (51%) were annotated 
(CanFam 2.0), it is possible that important genes are missed. 
Gene expression profiles of the two HS forms were compared with 
each other to identify genes that are specific for each manifesta- 
tion. 

When comparing VHS and STHS, 191 probes were signifi- 
cantly differentially expressed, listed in Additional Table $1. From 
the selection of most significant (P<0.003), unique probes with a 
threshold of log2 fold change of 1.5 (n= 19), eight were excluded 
for either comprising chromosomal regions (n = 3) or as clones that 
did not align (n=5). The eleven genes that remained, are 
visualized in a heatmap (Figure 1). Their potential involvement 
in tumor development or behavior was considered by a literature 
research. All microarray gene expression data were deposited in 
the public data repository GEO (accession number GSE45832). 

QPCR confirmed and quantified the differential expression of 
the 11 genes selected from the micro-array data. As a result, C6, 
S100A12, S100A8, PHYH and IGF were up-regulated and TATLI, 
CLEC12A, CD9, CCL5 and ASPN were down-regulated in VHS 
compared to STHS. Only for three gene products was the 
difference in expression significant: C6 (P=0.038), CLECI2A 
(P= 0.026) and CCL5 (P= 0.0069) (Figure 2). 

For technical reasons, no qPCR data could be obtained for 
CAMR2A. 

Our observations of CLECI]2A and CCL5 and those made in 
several human cancer types, make it conceivable that these gene 
products play a role in HS. In humans, C chemokine ligand 5 
(CCL5) functions as one of the natural ligands for the CC 
Chemokine Receptor 5 (CCR5). It mediates chemotactic activity in 
immune cells including monocytes and dendritic cells [36] CCL5 
[37,38] and CCR5 [39] promote breast cancer invasiveness and 
metastatic potential, while CCR5 inhibition abrogates this [39]. 
For inflammatory breast cancer, CCL5 is considered to constitute a 
prominent part of a poor prognosis signature [40]. Also in human 
colorectal carcinoma CCL5 appears to stimulate cancer progres- 
sion [41]. Furthermore, the CCL5/CCR5 axis has been shown to 
promote cell motility in human osteosarcoma [42]. Thus, for these 
human cancers CCL5 appears to be associated with a metastatic 
phenotype. Yet, in our study, CCL5 was expressed at lower level in 
VHS than in STHS. Considering the more rapid metastatic 
nature of VHS, this is surprising, and this finding needs 
confirmation (also at protein level) in future studies. Still, a 
decrease in CCL5 expression has been described in other human 
malignancies such as colon carcinoma when compared to normal 
tissue [43]. An alternative explanation for a reduction in CCL5 
expression in visceral as compared to localized HS could be that a 
reduction in expression of CCL5 could protect against immuno- 
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Table 2. QPCR primers for genes of interest (based on Microarray analyses) and reference genes (efficiencies varied between 95.0 


and 104.8%). 


Amplicon Annealing 


Gene name ENSID Primer sequence Forward/reverse size (bp) T (°C) 

CLEC12A ENSCAFG00000025113 AAATGCCAGCCTGTTGAC F 110 61 
TGGTAATCTCTGTCATACTTGGG R 

CCL5_CANFA ENSCAFG00000018171 TATGCCTCAGACACCACAC F 119 63 
GACAAAGACGACTGCTGG R 

ASPN ENSCAFG00000002307 CCACGAGTCAGAGAAATACAC F 135 59 
GGCAGAAGTCATTCACTCC R 

CD9 ENSCAFG00000015172 TGTGCTGTCATCCATCAC IF 118 57 
TGCCAAATATCATCACTACGG R 

TKTL1 ENSCAFG00000019451 ATGAGATACAAACAGGAGGAC IF 136 61 
CCAGTATATGCCATCCCAC R 

C ENSCAFG00000018598 CCGTTGTGATTGACTTTGAG IF 88 61 
CTTTCTGAGGTTGTTCCGT R 

S100A12 ENSCAFG00000023324 AAAGGGTGAGATGAAGCAG IF 156 61 
CACAACAGAAACCAGGGA R 

IGJ ENSCAFG00000002911 CCTTCTCCCGATGATCCT IF 120 63 
GGTACACAAATTTCGTTCTCAC R 

S100A8 ENSCAFG00000017557 GTTTACCACAAGTACTCCCTG F 148 63 
CCATCGCTATTGACATCCA R 

PHYH ENSCAFG00000023349 CTGAAGCCACACGATTATCC iF 112 58 
TCTCCTTTCTCCATCACGA R 

HPRT ENSCAFG00000018870 AGCTTGCTGGTGAAAAGGAC F 104 56 
TTATAGTCAAGGGCATATCC R 

RPS19 ENSCAFG00000028485 CCTTCCTCAAAAAGTCTGGG F 95 61 
GTTCTCATCGTAGGGAGCAAG R 

RPL8 ENSCAFG00000001677 CCATGAATCCTGTGGAGC IF 64 55 
GTAGAGGGTTTGCCGATG R 

SRPR ENSCAFG00000010474 GCTTCAGGATCTGGACTGC iF 81 61 
GTTCCCTTGGTAGCACTGG R 

RPL13 ENSCAFG00000019840 GCCGGAAGGTTGTAGTCGT iF 87 61 
GGAGGAAGGCCAGGTAATTC R 

GUSB ENSCAFG00000010193 AGACGCTTCCAAGTACCCC iF 103 62 
AGGTGTGGTGTAGAGGAGCAC R 

GAPDH ENSCAFG00000024323 TGTCCCCACCCCCAATGTATC F 100 58 
CTCCGATGCCTGCTTCACTACCTT R 

B2MG ENSCAFG00000013633 TCCTCATCCTCCTCGCT F 85 61 
TTICTCTGCTGGGTGTCG R 

RPS5 ENSCAFG00000002366 TCACTGGTGAGAACCCCCT F 141 62.5 
CCTGATTCACACGGCGTAG R 


Genes identified using microarray as being significantly different comparing Soft Tissue Histiocytic Sarcoma (STHS) and Visceral Histiocytic Sarcoma (VHS): C-type lectin 
domain family 12, member A (CLEC12A); C-C motif chemokine 5 Precursor (Small-inducible cytokine A5) (CCL5_CANFA); Asporin (ASPN); CD9 molecule (CD9);Transketolase-like 
1, transcript variant 1 (TKTL1); Complement component 6, transcript variant 1(C6); $100 calcium binding protein A12 (S100A12); Immunoglobulin J polypeptide (IGJ); S100 
calcium binding protein A8& (S100A8); Phytanoyl-CoA dioxygenase, peroxisomal like (PHYH) Reference genes primers for q PCR: Hypoxanthine phosphoribosyltransferase 

(HPRT), Ribosomal protein $19 (RPS19) ribosomalprotein L8 (RPL8), Signal recognition particle receptor (SRPR), Ribosomal protein L13, (RPL13), glucuronidase, beta (GUSB), 
Glyceraldehyde-3-Phosphate Dehydrogenase (GAPDH), Beta-2-Microglobulin (B2MG), 40S ribosomal protein $5 (RPS5), Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), 


B-2-microglobulin (B2MG), Ribosomal protein $5 (RPS5). 


doi:10.1371/journal.pone.0098258.t002 


surveillance [44] and hence be related to more aggressive behavior 


of HS. 


CLEC12A expression was found to be significantly lowered in 
VHS compared to STHS. CLEC/2A (or MICL) is considered a 
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negative regulator of granulocyte and monocyte function [45]. 


Activation of myeloid cells and recruitment to sites of inflamma- 


tion — but not increase or decrease in the level of differentiation — 
was accompanied by reduced expression [46]. Whereas normal 
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CAMK2A 


$100A12 


$100A8 


PHYH 


Figure 1. Microarray- based heatmap of 11 genes. C-type lectin domain family 12, member A (CLEC12A); C-C motif chemokine 5 Precursor (Small- 
inducible cytokine A5) (CCL5_CANFA); Asporin (ASPN); CD9 molecule (CD9);Transketolase-like 1, transcript variant 1 (TKTL1); Complement component 6, 
transcript variant 1(C6); $100 calcium binding protein A12 (S100A12); Immunoglobulin J polypeptide (IGJ); $100 calcium binding protein A8& (S100A8); 


Phytanoyl-CoA dioxygenase, peroxisomal like (PHYH). 
doi:10.1371/journal.pone.0098258.g001 


lymphocytes have no or low MICL expression, a relatively high 
expression of this gene in acute lymphoblastic leukemia was found 
to be associated with prolonged relapse-free survival [47]. How the 
reduced expression in VHS as compared to STHS relates to these 
findings remains to be determined, but a resulting increased 
migratory capacity as present in VHS could be an explanation. 

Complement component 6 (C'6) gene expression was significantly 
increased in VHS as compared to STHS. In acute leukemia’s an 
increase in circulating complement is common. As one early study 
has demonstrated increased expression of complement in mono- 
cytes by conditioned media of leukemic cells [48]. At present, a 
straight forward hypothesis on the functional consequences of the 
variation in C6 expression in the different forms of HS is not easily 
postulated but it might be associated with a reaction of the innate 
immune system to the neoplasm and not a direct effect of the 
neoplastic cell population. 

The variation in expression of C’6, CCL5 and CLEC124, all three 
members of the immune response, is worthy of follow up 
investigations, with focus upon their character as deranged 
histiocytes and may include comparisons with the canine reactive 
histiocytic diseases [1]. 

When comparing the results of qPCR and microarray the 
difference in the expression of some genes between STHS and 
VHS did not attain statistical significance. The variation between 
the two methods relates to the fact that a microarray experiment is 
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a semi-quantitative screening method and the qPCR quantitative. 
Many methodological factors can lead to a lack of correlation 
between array results and qPCR measurements [49]. The 
expression of CD9 (synonym MRP-1; motility related protein 1) 
in VHS was suppressed as compared to STHS at a level trending 
towards significance (P= 0.07) and further investigation of this 
gene in these sarcomas seems warranted in view of associations 
with reduced expression of this gene with tumor aggressiveness 
such as reported for other cancers. A lower expression of CD9 was 
found to be associated with the formation of bony metastases in 
studies using breast cancer cell lines [50,51]. Similarly, in human 
colon cancer patients, patients that lacked CD9 mRNA expression, 
had a worse prognosis than the cases that did express CD9 mRNA 
[52]. 

DAVID pathway analyses of the significantly differentially 
expressed genes in the micro array did not lead to the detection of 
altered expression of whole pathways. Also the genes chosen for 
qPCR confirmation did not appear to be related at the level of 
regulation. 

The alterations in gene function as detected in the current 
analyses, need follow up by subsequent investigations by use of 
antibodies - most still need to be developed and validated for use in 
the dog to examine an altered expression at protein level. 

When looking at tumor conditions in the human that share 
features with canine HS, none of the eleven genes for which 
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Figure 2. Quantitative PCR results. The upregulation or downregulation of selected genes in STHS (Soft Tissue Histiocytic Sarcoma) and VHS 
(Visceral Histiocytic Sarcoma). The thick black line represents the median (50th percentile) and also the first and third quartile (25th and 75th 
percentile respectively) are displayed.Three genes are significantly differentially expressed: C6 was up-regulated when comparing the more 
aggressive visceral histiocytic sarcoma to the soft tissue form, and CLEC12A and CCL5 were down-regulated when comparing the more aggressive 
visceral histiocytic sarcoma to the soft tissue histiocytic sarcoma. Abbreviations: C-type lectin domain family 12, member A (CLEC12A); C-C motif 
chemokine 5 Precursor (Small-inducible cytokine AS) (CCL5_CANFA); Asporin (ASPN); CD9 molecule (CD9);Transketolase-like 1, transcript variant 1 (TKTL1); 
Complement component 6, transcript variant 1(C6); $100 calcium binding protein A12 (S100A12); Immunoglobulin J polypeptide (IGJ); $100 calcium 
binding protein A8 (S100A8); Phytanoyl-CoA dioxygenase, peroxisomal like (PHYH). 

doi:10.1371/journal.pone.0098258.g002 


altered function was observed in the present study, have been development of HS, irrespective of form, this study provides the 
recognized as aberrantly expressed in micro-array/qPCR studies most comprehensive database to date of genetic variations in the 
in Langerhans cell histiocytosis [53]. two most common forms of HS, namely VHS and STHS. Using 

fold-change analysis, it reveals genetic variations not previously 
Conclusions associated with these two forms. 


os: : : ; On the basis of quantitative differences in expression of C’6 (up- 
As a valuable addition to our previous study, in which we were regulated in VHS versus STHS), CLEC/2A and CCL5 (down- 
able to provide evidence for involvement of several genes in the regulated in VHS versus STHS) were associated with each 
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subtype. Down-regulation of CLEC/2A in VHS, the more 
aggressive form of HS, is in line with previously published 
observations that such reduced function facilitates migratory 
capacity of myeloid cells in humans [46]. Down-regulation of 
CCLS is in line with several studies in human cancers [37,47], 
however contradictory to others [43]. Further investigations 
should focus on changes of gene function at protein level and a 
comparison of these histiocytic malignancies in dog and human. 


Supporting Information 


Table S1 All 191 probes significantly differentially 
expressed. Of four of the 191 probes, the annotation could 
not be traced back. Of the remaining probes, 28 have a ‘genomic’ 
location for which no gene could be mapped at this moment. 159 
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Probes have a gene-name, of which 142 are unique. From the 
selection of most significant (P<0.003), unique probes with a 
threshold of log2 fold change of 1.5 (n= 19), eight were excluded; 
either for comprising chromosomal regions (n= 3) or clones that 
did not align (n=5, marked with *). The eleven genes that 
remained, are written in bold and italic. 
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